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Abstract. Computer algebra methods are applied to investigation of spectral 
asymptotics of elliptic differential operators on curved manifolds with torsion and 
in the presence of a gauge field. In this paper we present complete expressions for 
the second coefficient (-E2) in the heat kernel expansion for nonminimal operator 
on manifolds with nonzero torsion. The expressions were computed for general case 
of manifolds of arbitrary dimension n and also for the most important for E2 case 
n — 2. The calculations have been carried out on PC with the help of a program 
written in C. 



1 Introduction 

Determination of the internal structure of an object via the spectra of differ- 
ent radiations and waves around the object is one of the archetypal problems 
in physics. More restricted mathematical version of this problem may be 
formulated as follows. A manifold (bundle) equipped with such structures as 
metric, curvature, torsion, gauge fields etc. and an elliptic (pseudo)differential 
operator acting on this manifold are given. What information about the man- 
ifold can one obtain studying the spectral properties of the operator? M. Kac 
phrases the problem in the evocative title of his paper "Can one hear the 
shape of a drum?" Q . The answer to this radical question is negative. In 
1964, J. Milnor Q found a pair of isospectral (with respect to the Laplace 
operator) but non- isometric tori in dimension sixteen, and recently, in 1999, 
D. Schiith ||^ constructed continuous isospectral families of metrics on the 
product of spheres S** x S*^ x S^l^ Nevertheless, many global geometric invari- 
ants, such as dimension, volume, and total scalar curvature, are known to be 
spectrally determined. Moreover, various manifolds such as round spheres of 
dimension less than or equal to six and 2-dimensional flat tori are uniquely 
determined by the spectra of the Laplacian acting on them. 

One of the most constructive approaches to study the spectral properties 
of operators on manifolds is investigation of the heat kernel expansion. This 
approach can be described briefly as follows. Starting with an elliptic operator 

^ After the Milnor's result many examples of multiply connected isospectral man- 
ifolds have been constructed, but the Schiith's construction is the first example 
of closed simply connected isospectral but non-isometric Riemannian manifolds. 
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A of the order 2r, acting on a bundle whose base is a compact close n- 
dimensional manifold M, and introducing an additional "time" variable t 
one can construct the heat operator A — Then one can compute the 
short-time asymptotic expansion of the diagonal elements of the kernel of 
this heat operator: 

(x|e-*^|x) ~ ^ E^(x|A)t — , t^+0. (1) 

m>0 

The coefficients EmixlA) in this expansion are spectral invariants of the op- 
erator A, and encode information about the asymptotic properties of the 
spectrum. These coefficients are called the heat invariants or heat kernel 
coefficients. They are also widely known under the names Hadamard co- 
efficients Hadamard-Minakshisundaram-DeWitt-Seeley, DeWitt-Seeley- 
Gilkey (HMDS or HAMIDEW, DWSG) coefhcients, according to papers of 
these authors ■ The heat invariants Em are of fundamental importance 

in quantum field theory, quantum gravity, spectral geometry and topology 
of manifolds. Many quantities of interest (such as the effective action. Green 
function, anomalies in quantum field theory P,|7|,|8|,p|jlc[| , the indices of elliptic 
operators and the invariants of manifolds in spectral geometry [pT|J^Jl],[l2[| ) are 
expressed in terms of the heat invariants. 

Most papers devoted to computation of the heat invariants deal with so- 
called minimal operators whose leading term is a power of the Laplacian 
and symbol is a scalar (w.r.t space-time indices). A typical example of such 
operators is 

A = -n+x. 

Here □ = g^'^D^D,^, is a covariant derivative including generally differ- 
ent connections (affine and spinor connections, gauge fields), AT is a matrix 
in internal space, i.e., an operator acting in sections of the bundle. For min- 
imal operators there are efficient enough methods for computing the heat 
invariants. 

In this paper we consider nonminimal operator of the form 

Af"" = -g^"' □ +aD^'D'' + X^"', (2) 

where X'^'^ is a tensor field (bundle indices are assumed implicitly), a is a 
scalar parameter which should satisfy to the condition a < 1 for the positive 
definiteness and, hence, for the ellipticity of operator (||). In recent years 
special cases of operator (H) have been encountered by physicists studying 
the quantization of gauge and gravitational fields in arbitrary gauges ||l0| , [r3| . 
For example, the quantization of Yang-Mills field in an arbitrary covariant 
background gauge leads to the operator 

AH - -5^, n^'-i^^-l^ Dl'^Df - 2r^Gl,, 

^ It was Hadamard who introduced these coefficients for scalar operator A already 
in 1923 and established their essential properties. 
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where Z?^ is a covariant derivative containing the external field potential , 
Gfj_i, is a corresponding field strength, /'^'"^ are the structure constants of a 
corresponding Lie algebra and a is a scalar (gauge) parameter. Another ex- 
ample: the quantization of electro-magnetic field in an external gravitational 
field leads to the operator [ p^[jT5[ | 

A,,^ = -g^^ □ - - 1^ Df,D^ + Rf,^ 

(for an analogous operator in quantum gravity see fl^ ). 

The torsion is defined as antisymmetric part of affine (or linear) connec- 
tion 

where I^'*'^^ are connection coefficients. The Einstein's General Relativity (see 
e.g. jl^) is based on a special connection called Levi-Civita connection, i.e., 
symmetric and compatible with metric affine connection. This connection can 
be expressed completely in terms of metric and is torsionless. The General 
Relativity well describes the interaction of the matter with the gravity as far 
as macroscopic bulk matter is considered. However, on the microscopic level, 
where the elementary particles posses such quantum property as spin, it seems 
necessary to take into account the infiuence of spin on the geometry of space- 
time. To describe the interaction of spinning particles with the gravitation, 
a gravitation theory should include the non-vanishing torsion. In 1922 Elie 
Cartan first pointed out |l^Jl^ that there is no a priori reason to assume 
an affine connection to be symmetric in the context of General Relativity. 
He proposed also a theory of gravitation with torsion which development is 
known now as the Einstein-Cartan theory. The torsion arises naturally in 
the different (based on Poincare and affine groups) gauge theories of gravity 
developed in the recent years (see Refs. pj| , pO[ ). Moreover, all kinds of modern 
superstring theories (for the recent review see e. g. [^), which allow to 
deduce the properties of space-time, also predict, along with the metric, the 
existence of torsion. 

In psf we computed the heat invariants for operator (^ up to E4 but 
for manifolds without torsion. In this paper we consider more general (and 
computationally much more difficult) case of manifolds with torsion. 



2 Algorithm and Implementation 

The algorithm we use was developed by V. Gusynin |^,^. This algorithm is 
based on the covariant generalization of the pseudodifferential calculus given 



by Widom |26 . The main advantage of this algorithm is its universality. It 
can be applied to the wide class of pseudodifferential operators, in particular, 
to the nonminimal and higher-order operators intractable by other methods 
such as the DeWitt ansatz Iql for heat kernel matrix elements. 
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The algorithm has the following main features. For a positive elliptic 
operator A the spectrum of which lies inside a contour C, the heat operator 
exp(— iA) can be expressed in terms of the resolvent {A~X)^^ via the formula 



Jc 27r 



The pseudodifferential calculus method uses the following representation for 
the matrix elements of the resolvent 

A — A J [2TTy''y'g(x') 

where cr{x, x' , fc; A) is an amplitude, l{x, x' , k) is a (real) phase function which 
is a biscalar with respect to general coordinate transformations, k is a wave 
vector. 

The resolvent satisfies the equation {A — X)G = 1 which leads to the 
equation for the amplitude: 

{A{x, Df, + iD^l) - A)cr(x, x',k; A) = I{x, x'), (5) 

where I{x, x') is a transport function having both bundle and Lorentz indices. 

In the pseudodifferential calculus, it is assumed that in the flat space the 
phase function has the form I = {x — x')^k^^ . The covariant analogue of the 
linearity of the function I is based on the requirement that all higher-order 
symmetrized covariant derivatives of l(x^ x\ k) vanish at the points x — x' , i. 
e., satisfy the infinite set of relations ||2q| : 

[{D^,...D^Jl] = Q, m>l, (6) 

where {. . .} means symmetrizing in all indices, and [. . .] means transition 
to coincidence limit [x = x'). In an analogous way, the covariant transport 
function should satisfy the relations: 

[{D^,...D^„^}I]^0, m>l. (7) 

Equations (|^) and together with the "initial conditions" [I] = 0, [D^^l] = k^ 
and [/] = l (unit operator) allow one to compute the coincidence limits 
for nonsymmetrized covariant derivatives [I?^^ . . . -Dp„J] and [Dn-^ . . . Dn^I]. 
These nonsymmetrized derivatives are obtained directly from (ph and (R) by 
reducing all terms to a unified index ordering with the help of the Ricci iden- 
tity. The resulting expressions are universal polynomials in the torsion T^^jy, 
curvature tensor R^f_i,^ri, gauge curvature W^i, and their covariant derivatives. 
In fact, once computed and stored the coincidence limits [-D^^ . . . D^j^^l] and 
[iJ^j . . .D^^I] can be used in many calculations for different operators A. 
The functions l{x,x',k) and I{x,x'), introduced with the help of formulas 
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and play an important role in the covariant pseudodifferential cal- 
culus called also intrinsic symbolic calculus p6t . In fact, just these universal 
functions manifest the geometric properties of a base manifold and a bundle. 

Expanding the amplitude a in degrees of homogeneity of k: 

00 

m—l 

we obtain the recursion equations for am from equation (^). For example, for 
operator (||) these recursion expressions take the form 

^^^CTiA^ + i [-g^^iO I + 2D'^lDrj) + a{D''D^l + D>'ID^ + dHd^")] 
X (ToA^ = 0, 



Af'^drnXu + i [-.g''^(n I + 2DilD,j) + aiD^'DH + D^HD^ + DHD^")] 
X <J(m-i)Xu + (-.9^^ □ +aD^D^ + X^^)a^m-2)Xu =0, m > 2, 
where the matrix 

^ gi^^i^D'^iD^^l - A) - aD'^lD'^l 

is the principal symbol for operator (H). Solving the recursion equations we 
obtain am- The heat invariants are expressed in terms of integrals of the 
coincidence limits [am] ■ 

/i-inh. r ;j\ 

(2^0V^ Jc ^^"'[^"K^,fc, A) ^ J[[am]). (8) 

The integrals in (^) can be expressed in terms of gamma and Gauss hyper- 
geometric functions for a wide class of operators A. The typical integral of 
terms of the coincidence limit [am] takes the form 

t/ k'^^k^-^ . . . k^2B \ _ 

V(fc2'-- A)'[(l-a)fc2'-- A]™y ~ 

r{{p + s + n/2)/r) ^, , , N 

where g{pi...p2s} ^ symmetrized sum of products of metric tensors. Using 
the fact that to and I are whole numbers, one can express the hypergeometric 
function in in terms of elementary functions with the help of the Gauss 
relation 

a(l — z)F(a-(- 1, 6; c; z) — [c— a)F{a— 1, 5; c; z) + {2a — c — az + bz)F{a, 6; c; z), 

(9) 

^ The existence of these functions has been proved in |26| . 
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and using then the formula |2 

■(l-fo)m-l 



fc! 

fc=0 



TO=1,2,..., m-6^1,2, 



where (a)^ = a{a + 1) . . . (a + fc — 1) is the Pochhammer symbol (shifted fac- 
torial). During simplification of tensor expressions we use various symmetry 
properties of the tensors -R^,,;^!/, T^^i^, W^^, and also the Ricci identity 

1=1 

k 

/ ^>A"''/^Ai...A,_iQA,+i...Afc + -'- A"'-^a'f'Ai...Afc + '^'^A'f'/^Ai...Afc' 

i=l 

the Bianchi identities for both affine and gauge curvatures 

+ W^aAT^0^ + T^^aT^^ + W^xT'^c.p = 0, 
and the cyclic identity 

R'^-^S + R'^513 + ^5/37 + DpT^-^s + DjT°'sp + DsT°'pj 
+ PX-L 7(5 + J 7AJ + J i5aJ Pf — U. 

The above algorithm has been implemented in the C language. The C code 
of total length about 11000 lines contains about 250 functions for different 
manipulations with tensors and scalars. These functions are gathered into 
two programs DWSGCOEF and COLIM. 

The COLIM program computes coincidence limits of the l{x,x',k) and 
I{x, x') functions and writes them to the disk. Once computed and stored^ the 
coincidence limits, being universal functions, can be used in many calculations 
for different operators A. 

The DWSGCOEF program computes Em coefficients by the following 
steps: 

1. Reading input information (operator, order m, etc.) 

2. Computing a set of asymptotic operators for constructing recursion equa- 
tions. 



* For the operators of different tensor ranks A,A'^'', . . . the coincidence limits for 
the functions /, J^", . . . should be computed separately. 
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3. Computing Um with the help 0/ the recursion equations. 

4. Taking the coincidence limit [um] ■ 

5. Integrating [am] to obtain the coefficient Em- 

6. Substituting tensor expressions for [D^^ . . . D^^l] and [D^-^ . . . D^^I] into 

Em- 

7. Reducing hypergeometric to elementary functions in the scalar coefficients 
( Ci in the formulas of Section |^ including in the heat invariants in 
the case of nonminimal or higher-order operator. Eliminating possible 
linear dependencies among these scalar coefficient^ to make the resulting 
formulas as compact as possible. 

8. Output Em (and its Lorentz trace in the nonminimal case). 

To cut down the swelling of the intermediate expressions, we use term-by- 
term strategy, i.e., the most cumbersome Steps are applied consecutively 
to single terms of am generated during the execution of Step ^. 



3 Heat Invariant E2 

We present here the full expression for the coefficient E2 and also its trace 
with respect to Lorentz indices for nonminimal operator (Q) on a curved man- 
ifold with the torsion and gauge field manifested itself in the gauge curvature 
W^^. We consider the case of arbitrary dimension n and also the most im- 
portant^ for E2 case n = 2. In the below formulas the indices a, /3, 7 and /i, i/ 
are dummy and free, correspondingly. We use the following definition for the 
torsion trace: = T"ct^. 



3.1 Full Expression 



E2 = (47r)-- ^-CiX"" - C2 {X"" + g'^'X^") + C3 i^W"" -t- -Da^T""" 
19 \ 



Usually there are many dependencies among the scalar coefficients which are not 
seen in terms of hypergeometric functions, i.e., quite different hypergeometric 
expressions may be reduced sometimes to the same elementary function. 
The Atiyah-Singer index of an elliptic operator on a manifold of dimension n 
can be expressed in terms of an integral of E„ over the manifold. 
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Coefficients Cj in arbitrary dimension n: 

Ci = — ^— I (1-a)"^ (San - 6a + 4n + 4) + an^ - 2an^ 

a{n — 2)n(n + 2) 

-San + 6a - 4n - 4} , 
= ^^^^-^1^^ { (1 -a)- 1 (an + 2a - 4) + an - 2a + 4} , 

C3 = ^_J_^|(i_a)i-t(an-8) + 3an-8a + 8|, 
a{n — 2)n ' ^ ' > 

Ci = TT-, ^ l(l-a)"^ (na^n^ + 34a^n - 17an^ - 168an 

6a(n- 2)n(n + 2) "-^ ' ^ 

-268a + 140n + 256) - SSan^ + 40on + 268o - 140n - 256} , 

C5 = — ^ l(l-a)"-^ ilba^n^ + ?,Qa^n - Iban^ - 152on 

6a(n-2)n(n + 2) ' ^ 

-244a + 116n + 256) - 43an^ + 24an + 244a - 116n - 256} , 

Ce = r l(l-a)"* (a^n^+2a^n-aV-20a^n-36a^ 

6a2(n-2)n(n + 2) "-^ ' ^ 

+12an + 96a - 48) + a^n^ - l&a^n + 36a^ + 12an - 96a + 48} , 

C^ = -4- 77 -|(l-a)-7 (-aV - 6aV + o^n^- 

6a2(n-2)n(n + 2)(n + 4) ' ^ 

8a3n + SOa^n^ + 152a2n-24an2 + 192a2-240an-576a+ 144n 
+288)- Ta^n^ + 18a2n^ + 64a^n-48an2-192a^ + 96an + 576a 
-144n-288}, 

= e,2(,_,'),(,^,) {(l-«)^-^ {a^n^ + 2a^- - 24an - 48a+ 144) 

-7a^n^ + SAa^n - 48a^ - 48an + 192a - 144} , 

Cg = --^ -7 |(l-a)"5 (a^n^ + 6a^n + 8a^ - 12an 

a^{n-2)n{n + 2){n + 4)^^ ' ^ 

-48a + 48) - a^n^ + 6a^n - Sa"^ - 12an + 48a - 48} , 

Cio = TTT^TT rr4 ?:t7 TT (-a^n^ - 6aV + a^n^- 

12a2(n-2)n(n + 2)(n + 4) "-^ ^ ^ 

8a^n + SOa^n^ + 152a^n + 1920^ - 192an - 768a + 576) - a^n^ 
-6a^n^ + 88a2n - 192a^ - 96an + 768a - 576} , 

Cn = —r- i-- — |(l-a)^"2 (-9aV - 18a^n + 76an + 152a 

3a"^(n — 2)n(n + 2) ^ 



-24) + 2 (-ISa^n^ + 60^71 + 760^ - 32an - 88a + 12) } , 

—r- ^— — {2(l-a)i-t (-5a^n^ - lOa^n + 38an 

3a2(n-2)n(n + 2) ^ ^ ^ ^ 

+12) -Sla'^n'^ + 26a2n + 1520^ - 88an - 128a - 24} , 



Heat Invariant E2 



1 



n 



6a(n- 2)n(n + 2)(n + 4) ^ ^ 

- 18a^n^ - 80a^n + 12an^ - 96a^ + 72on + 96n - 48n + 96) + a^n^ 
-ISa^n^ + Sa^n + I2av? + 960^ - 120an - 96a + 48n - 96} , 

= 2a^(n-2)n(n + 2) ' ''^'^ + + " 

-a^n^ - 2a^n + IGa^ - 32a + 16} , 

Ci5 = T-, (-aV-2a^n+an^+8an+12a-24) 

6a(n — 2)n[n + 2) ^ ^ ' 

+an^ - an^ - 12a + 24} , 

C16 = TTT^TT ^^4 7T {(l-^)"'^ (-aV - 6a3n2 + a^n^ 

12a2(n-2)n(n + 2)(n + 4) ' ^ 

-8a^n + Ga^n^ + Sa'^n + A8an + 192a - 288) + a^n^ + 3a^n^ 
+20^71^ - ASa^n + 96an - 192a + 288} , 

Ci7 = T-r^TT ^rr^ jt (-a^n^ - 6a^n^ + a^n^ 

24a2(n-2)n(n + 2)(n + 4) "-^ ^ ^ 

-8a^n + SOa^n^ + 152a2n + 192a^ - 192an- 768a+ 576) +a^n^ 

+3a^n^ - lOa^n^ + 72a2n - 192a2 - 96an + 768a - 576} . 

Coefficients Cj in the dimension n = 2; 

31n(l-a) 2-a 



4a 8(l-a)' 



^^_ ln(l-a) ^ 2-a 



a 2(l-a)' 
(a-4)ln(l-a) , ^ 

= 2-a + 

_ (17a -67) In(l-a) 137a - 134 
^ ~ 12^^ ^ 24(l-a) ' 

(15a-61)ln(l-a) 119a- 122 
12^^ ^ 24(l-a) ' 

^ (a^ - 9a + 6)ln(l-a) (3a-2)(a-2) 



5 



Ct 



12a2 8a(l-a) ' 

(a- - 12r/ + 12)lii(l-a) 2a'^ - 9a + 6 
12^? 6a(l-a) ' 



_ (a^ - 12a+ 18)ln(l-a) a- 6 
12^:5 

(a - 2)lii(l-r/) «--12r/ + 12 
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Cio 

Cl2 
Cl3 
Ci4 

Cir 



_(a2 - 12a+ 12)ln(l-a) _2a^ -9a + 6 
24a? 12a(l-a) ' 

(9a^ - 38a + 3)ln(l-a) _ 73a -6 
6a2 12a ' 

(10o2-38a-3)ln(l-a) 79a + 6 
6^ 12a ' 

(a-6)ln(l-a) 2a - 3 

12^^ ^ 6(l-a)' 

(o^ -4o + 2)lii(l-a) 3a -2 

4a~' 

_(a- 3)ln(l-a) _ 7a - 10 

12a 24(l-a)' 
(g^ -6)lii(l-a) 3a^ + a - 6 

24^2 24a(l-a) ' 

(o^ - 12a+ 12)ln(l-o) Sa^ - 10a + 6 



48a2 24o(l-a) 
3.2 Lorentzian Trace 

tvi^E^ = (47r)-t {-CiX„" + C2R - C^T^p^T'^"'' + dTaff^T^''^ 

+C5DoT« + (C4 + C5)T„T"} 

Cj for the trace in arbitrary dimension n: 

(1 — a)~^+n — 1 {l — a)~^{—an + n + 6)+n^ — n — 6 

= , (^2 = ^ , 

n on 

C3 = — — — |_(i_a)-f (a^n^ + 2a^n - an^ - 26an - 48a + 96) 



24on(n 


+ 2) 


+an^ 4 


- an^ - 


1 




12an{n 


+ 2) 


-an^ - 


an^ - 


1 





C4 = — — — {(1-a)--^ {a^n^ + 2a?n - an^ - Uan - 24a + 48) 



C5 = — — |-(l-a)"^ ia^'n? + Aa^n - an^ - lOan - 36a + 48) 

6a(n — 2)n 

+an^ - 3an^ + 14an - 36a + 48} . 
Ci for the trace in the dimension n = 2: 



2-« ^ _ 2+a ^ _ 1 ^ _ 1 
2(1^' - 6(1^' - 12' - ~6' 

(a-4)ln(l-a) 11a- 14 



2a 6(1 -a) ' 
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4 Conclusion 

The program computes E2 with torsion for operator (||) rather easily (about 
10 sec on a Pentium-75 PC). Unfortunately, computational complexity of 
the problem under consideration is very high. For example, the timings for 
torsionless computations of E2 and E4 for the same operator (||) are < 1 sec 
and 4 h 5 min, correspondingly. It is clear that the inclusion of the torsion 
increases the computational efforts considerably and the computation of i?4 
with torsion may take too much time. Another problem is the volume of the 
resulting expressions. There are two ways to handle this problem. First of 
all, some work is needed for developing of algorithms for further reduction of 
large tensor expressions. However, due to the natural complexity of the heat 
invariants, one can not hope to make the higher-order invariants tractable by 
hand. Thus, the methods for automatic usage of these invariants should be 
elaborated. 
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